Negative thermal expansion in NbF_3 and NbOF_2:A comparative theoretical study

Thermal expansion control is always an obstructive factor and challenging in high precision engineering field. Here,the negative thermal expansion of Nb F_3 and Nb OF_2 was predicted by first-principles calculation with density functional theory and the quasi-harmonic approximation(QHA). We studied the total charge density, thermal vibration, and lattice dynamic to investigate the thermal expansion mechanism. We found that the presence of O induced the relatively strong covalent bond in Nb OF_2, thus weakening the transverse vibration of F and O in Nb OF_2, compared with the case of Nb F_3.In this study, we proposed a way to tailor negative thermal expansion of metal fluorides by introducing the oxygen atoms.The present work not only predicts two NTE compounds, but also provides an insight on thermal expansion control by designing chemical bond type.     

NUMERICAL SIMULATION OF FIBROBLAST DUROTAXIS MIGRATION INDUCED BY STIFFNESS GRADIENT OF SUBSTRATE1)

建立了一种细胞趋硬性迁移的理论模型和有限元分析框架,为连续变刚度人工基质的试验设计提供理论依据。考虑了细胞体的黏弹性属性,以及细胞与基质间的配受体动态反应过程,并以配受体合成时间为时间步长,将细胞运动方程化为静力学形式进行求解。对有限元过程提出一种动约束,便于消除其结构矩阵的奇异性。结果表明,模型能够模拟黏着斑内部力的快速波动现象,细胞的运动速度与观测数据一致,可有效模拟20,h以上的长时程问题。     

COMPUTER-ASSISTED TEACHING REFORM FOR THEORETICAL MECHANICS 1)

本文依照四川大学某理论力学教学班实际教学改革的经验,通过具体举例说明计算机技术在理论力学教学中尚未开发的巨大潜力,减小数学运算使学生重点关注理论力学基本原理,批量求解使学生提前接触实际应用案例,普适性求解使学生深化对处理多变量问题的认识。同时本文还探讨了理论力学教学改革的方向和目的,即把教学重点放在普适的教学方法上,借助计算机求解多刚体、多变量问题,以期有助于理论力学教材更好地适应时代发展。     

水化硅酸钙纳米压痕分子动力学模拟方法优化

针对水化硅酸钙纳米压痕模型忽视了压头与基底之间相互作用的问题，由尺寸差异引起的金刚石压头难以计算的问题，以及Wittmann模型无法得到实际接触面积的问题，提出了新的模型与计算方法．结合分子动力学方法，采用金刚石压头-Wittmann模型基底的组合方式构建无定形态水化硅酸钙纳米压痕试验模型．在建模阶段，考虑到压头模型与基底模型粒子间尺寸差异，提出了等比例替换模型，通过公式推导并就不同尺寸模拟结果验证了等比例替换模型的可行性．在计算阶段，提出了局部前处理的弛豫方法进行模拟．确定最大荷载位置处的接触面积为546 nm2，进而求出水化硅酸钙模型硬度H为0.84 GPa、折合模量Er为30.52 GPa．并通过纳米压痕试验，验证了模拟结果的准确性，证明了模型的科学性，对今后水化硅酸钙(C-S-H)纳米层面的模拟具有重要借鉴意义．     

Bridging the Green Gap: Metal–Organic Framework Heteromultilayers Assembled from Porphyrinic Linkers Identified by Using Computational Screening

In organic photovoltaics, porphyrins (PPs) are among the most promising compounds owing to their large absorption cross-section, wide spectral range, and stability. Nevertheless, a precise adjustment of absorption band positions to reach a full coverage of the so-called green gap has not been achieved yet. We demonstrate that a tuning of the PP Q- and Soret bands can be carried out by using a computational approach for which substitution patterns are optimized in silico. The most promising candidate structures were then synthesized. The experimental UV/Vis data for the solvated compounds were in excellent agreement with the theoretical predictions. By attaching further functionalities, which allow the use of PP chromophores as linkers for the assembly of metal-organic frameworks (MOFs), we were able to exploit packing effects resulting in pronounced redshifts, which allowed further optimization of the photophysical properties of PP assemblies. Finally, we use a layer-by-layer method to assemble the PP linkers into surface-mounted MOFs (SURMOFs), thus obtaining high optical quality, homogeneous and crystalline multilayer films. Experimental results are in full accord with the calculations, demonstrating the huge potential of computational screening methods in tailoring MOF and SURMOF photophysical properties.     

Ab initio and Monte Carlo studies of phase transitions and magnetic properties of YCrO3: Heisenberg model

By using the density functional theory (DFT) and Monte Carlo simulations (MCS) with the Heisenberg model, we have studied magnetic properties of the bulk perovskite YCrO₃. The exchange couplings of the Heisenberg model and the magnetic anisotropy are investigated. The 110 direction in the crystalline structure of the compound has shown the minimum energy, it is the easy magnetic direction. Using Monte Carlo simulations, the magnetizations behavior, the effects of system parameters and the critical exponents of the compound YCrO₃ are implemented. It is shown that the bulk perovskite YCrO₃ belongs to the 3D Heisenberg universality class.     

A Theoretical study on the charge and discharge states of Na-ion battery cathode material,Na1+xFePO4F

Na₂FePO₄F is a promising cathode material for a Na-ion battery because of its high electronic capacity and good cycle performance. In this work, first principle calculations combined with cluster expansion and the Monte Carlo method have been applied to analyze the charge and discharge processes of Na₂FePO₄F by examining the voltage curve and the phase diagram. As a result of the density functional theory calculation and experimental verification with structural analysis, we found that the most stable structure of Na_1.5FePO₄F has the P2₁/b11 space group, which has not been reported to date. The estimated voltage curve has two clear plateaus caused by the two-phase structure composed of P2₁/b11 Na_1.5FePO₄F and Pbcn Na₂FePO₄F or Na₁FePO₄F and separated along the c-axis direction. The phase diagram shows the stability of the phase-separated structure. Considering that Na₂FePO₄F has diffusion paths in the a- and c-axis directions, Na₂FePO₄F has both innerphase and interphase diffusion paths. We suggest that the stable two-phase structure and the diffusion paths to both the innerphase and interphases are a key for the very clear plateau. We challenge to simulate a nonequilibrium state at high rate discharge with high temperature by introducing a coordinate-dependent chemical potential. The simulation shows agreement with the experimental discharge curve on the disappearance of the two plateaus. © 2018 Wiley Periodicals, Inc.     

Critical Assessment of the DFT + U Approach for the Prediction of Vanadium Dioxide Properties

In a previous study (Stahl and Bredow, Chem. Phys. Lett. 2018, 695, 28–33), we have studied structural, energetic, and electronic properties of two vanadium dioxide VO₂ polymorphs with modified global and range-separated hybrid functionals. Since hybrid methods are computationally demanding, we evaluate the computationally more efficient DFT + U method in the present study. We assessed the widely used Dudarev PBE + U approach with a literature value of the effective Hubbard parameter U_eff = 3.4 eV. We compared the PBE + U results for the two VO₂ polymorphs with our previous results, a self-consistent hybrid functional sc-PBE0, and the meta-GGA functional SCAN. It was found that the PBE + U method yields a strongly distorted monoclinic phase and does not reproduce the metal-to-insulator transition of VO₂ correctly, even with modified values of U_eff. On the other hand, sc-PBE0 and SCAN describe the relative stability and the electronic structure of both polymorphs correctly and also provide reasonable lattice parameters. The functional SCAN yields the optimal balance between computational efficiency and accuracy. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.     

双原子分子激发态势能的自旋相关局域Hartree-Fock密度泛函理论方法

基于自旋相关局域Hartree-Fock (SLHF)势函数，本文提出了一种计算双原子分子激发态势能的密度泛函理论(DFT)方法，并将该方法应用于和的激发态势能曲线的计算。在只考虑交换能的情况下，本文的DFT计算结果与文献中精确方法和Hartree-Fock (HF)方法的结果符合的非常好，说明采用SLHF势函数作为交换势的DFT方法是一个很好的计算激发态势能的方法。本文还计算和探讨了电子的关联势函数和关联能，发现传统的近似方法在较大核间距的情况下大大低估了电子的关联能.     

Identifying the viscoelastic properties of soft matter from the indentation response of a hard film-soft substrate system

The indentation technique is widely used in measuring the mechanical properties of soft matter at the microscale or nanoscale,but still faces challenges by these unique properties as well as the consequent strong surface adhesion, including the strong nonlinear effect, unclear judgment of the contact point, difficulties in estimating the contact area, and the risk of the indenter piercing the sample. Here we propose a two-step method to solve these problems: lay a hard film on a soft matter, and obtain the viscoelastic properties of this soft matter through the indentation response of this composite structure. We first establish a theoretical indentation model of the hard film-soft substrate system based on the theory of plates, elastic-viscoelastic correspondence principle and Boltzmann superposition principle. To verify the correctness of this method, we measure the mechanical properties of the methyl vinyl silicone rubber(MVSR) covered by a Cu nanofilm. Finally, we test the effectiveness and error sensitivity of this method with the finite element method(FEM). The results show that our method can accurately measure the mechanical properties of soft matter, while effectively circumventing the problems of the traditional indentation technique.